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Abstract 

We study on the spectrum and short-distance two-body force of holographic baryons by the matrix 
model, which is derived from Sakai-Sugimoto model in D0-D4 background (D0-D4/D8 system). The 
matrix model is derived by using the standard technique in string theory and it can describe multi- 
baryon system. We re-derive the action of the matrix model from open string theory on the wrapped 
baryon vertex, which is embedded in the DO- D4/D8 system. The matrix model offers a more systematic 
approach to the dynamics of the baryons at short distances. In our system, we find that the matrix model 
describe stable baryonic states only if C = < 2, where Uq^ is related to the number density 

of smeared DO-branes. This result in our paper is exactly the same as some previous presented results 
studied on this system as |27| . We also compute the baryon spectrum (fc = 1 case) and short-distance 
two-body force of baryons (k = 2 case). The baryon spectrum is modified and could be able to fit the 
experimental data if we choose suitable value for C,. And the short-distance two-body force of baryons 
is also modified by the appearance of smeared DO-branes from the original Sakai-Sugimoto model. If 
C > 2, we find that the baryon spectrum would be totally complex and an attractive force will appear in 
the short-distance interaction of baryons, which may consistently correspond to the existence of unstable 
baryonic states. 
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1 Introduction 


QCD as the underlying fundamental theory to physicists, has achieved great successes in particle physics and 
nuclear physics, however we still can not describe or predict the behavior of baryons or nucleons exactly. In 
high energy physics, it is well-known that nuclear physics remains one of the most difficult and intriguing 
branches. The key problem is that we have not yet understood much about the strong-coupling QCD 
since to study on the strong-coupling QCD is hopeless by using the techniques from perturbative quantum 
field theory directly. However about two decades ago, the discovery of AdS/CFT and gauge/string duality 
m 11 s] became one of the turning points. The AdS/CFT correspondence has been recognized as a 
promising framework to understand non-perturbative aspects of gauge field theory, thus it may be able to 
provide a new way to study on nuclear physics. 

There have been many applications or models of AdS/CFT for studying on strong-coupling QCD such as 
[3 H13 H H. Several top-down models for QCD inspired by AdS/CFT have also been proposed in recent 
years. The most well-known D4-D8/D8 brane system based on Witten’s work (TO] was built up by Sakai and 
Sugimoto [111I12| (Sakai-Sugimoto model). The Sakai-Sugimoto model consists of W D4-branes and Nf D8- 
branes as probes. The D4-branes are compactified on a cycle, offering color numbers, which represent pure 
Yang-Mills in the low energy effective theory. Supersymmetry is broken down by introducing the anti-periodic 
boundary conditions for fermions on the compactified cycle. Since only the low energy theory is concerned 
in this model, so the dynamics of background geometry produced by Nc D4-branes, can be described by 
Type II A supergravity on the large-Yc limit. The Nf species of massless flavored quarks are introduced by 
the embedding Nf pairs of probe D8/D8-branes. The flavor D8/D8-branes are connected at the IR region 
of the D4 solitonic solution, which corresponds to the geometrically broken chiral symmetry in the confined 
phase for the dual held theory. In this geometry, the low energy effective theory of light meson sector comes 
from the worldvolume theory on the connected D8/D8-branes. There is an other solution for the background 
geometry for this model, i.e. the D4 black brane solution, which corresponds to the deconhned phase for the 
dual held theory [13] (However, in fact it is less clear that what phase the deconhned geometry corresponds 
to in the dual held theory M)- As it is known that baryons are D-branes wrapped on non-trivial cycles 
[i5|[i6i[n], so in this model, baryons are identihed as D4-branes wrapped on a four-cycle which is named 
as the “baryon vertex”. And it has turned out baryons can be treated as a small instanton conhguration in 
the worldvolume gauge theory on the probe D8-branes [T5]. According to these viewpoints, there are many 
researches on baryons or nuclear matters by holography, such as the phase structure [Udiiiiniiii] and the 
interaction [23[531[21|- However some results from the Sakai-Sugimoto model is still not realistic to QCD. 
One of the most likely reasons may be that the Sakai-Sugimoto model is a theory with large W limit, but 
the real QCD is not. Therefore some generalizations or modifications of this model have been proposed such 
as and the backreaction of the flavor branes has also been considered recently as [5^. In our paper, 
we follow [25 and use the gauge/gravity duality to study the dynamics of multi-baryons from D0-D4/D8 
system by matrix model proposed in |24| . By using the matrix model, we calculate some basic properties of 
holographic baryons: the spectrum of baryons and the effective two-body potential both from Sakai-Sugimoto 
in D0-D4 background as an extension of 113 [I 3 . 

The original matrix model is the effective theory for DO-branes |1S| , which can also be understood as a dual 
description of D=II supergravity. As a generalization, the matrix model in [24], is proposed as the effective 
theory for baryons or nucleons by holography. In this matrix model, the rank of the matrix represents 
the number of baryons and k baryon branes produce U{k) symmetry for /c-body baryons. The diagonal 
elements of matrices represent the positions of k baryons while the off-diagonal elements are integrated out. 
Besides, the size of baryons are related to the classical values of a pair of complex k x Nf rectangular 
matrices, which describe the dynamics of the strings connecting the flavor branes and the baryon vertices 
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in low energy effective theory. With all together, it comes to the well-known Atiyah-Drinfeld-Hitchin-Manin 
(ADHM) matrix of instantons. So in our paper, we re-derive the matrix model from Sakai-Sugimoto in 
D0-D4 background by using standard technique in string theory. Such background geometry is produced by 
Nc D4-branes with Nq smeared DO-branes. In this background, we follow the original idea in USE!, and 
recall that baryons can be identified as D4-branes wrapped on the 4-cycle. Here we start to use D4’-branes to 
distinguish such a baryon vertex from those D4-branes who are responsible for the background geometry. In 
the large Nc limit, the dynamics of the open strings with both ends on the D4’-branes and which connects the 
D4’-branes and the D8-branes are also relevant, thus the dynamics of k baryons are described by U{k) gauge 
theory. And the theory on the D4’-branes is reduced to a 0+1 dimensional matrix model by considering only 
the zero modes along the on which the D4’- and the D8-branes are wrapped. Then it is clear that the 
matrix model is just the low energy effective theory for the baryon vertices. 

Our motivation for this paper is to study the holographic baryons by the matrix model derived from 
Sakai-Sugimoto model in D0-D4 background (i.e. D0-D4/D8 system). In this paper, since the appearance 
of smeared DO-branes could be able to modify the results about baryons from the original Sakai-Sugimoto 
model, thus our results may be more close to the realistic physics. On the other hand, accommodating many 
bodies of baryons by such a matrix model would be very easy, so we can even use this model to describe 
the interaction among multi-baryons. The outline of this paper is as follows. In Section II, we give a brief 
review of Sakai-Sugimoto model in D0-D4 background (D0-D4/D8 system). In Section III, we start to derive 
the matrix model in the D0-D4 background by considering D4’-D8 gauge theory compactified on S'^. Since 
we keep k baryons close to each other, so only the short-distance effects are relevant. It turns out that our 
matrix model can describe stable baryonic states only for ^ = Uq^/U^j^ < 2, where Uq^ is related to the 
number density of smeared DO-branes. And the constraint for stable baryonic states (C < 2) is exactly the 
same as m- And as two simple examples, we calculate the energy functions of static configurations with 
k = 1 and k = 2, for one and two ffavor(s). In Section IV, we use our matrix model derived from D0-D4/D8 
system to study baryon spectrum. We determine the size of the holographic baryons with the case of A: = 1, 
also for one and two flavor(s). Again it turns out that the baryon spectrum does make sense only for C < 2, 
otherwise baryons may not exist or be unstable. In Section V, we study the case of fc = 2 and use the 
instantons of ADHM matrix as data to calculate a baryon-baryon potential at short distance for two-flavor 
case. By integrating out the auxiliary gauge potential in 0+1 dimension, it also turns out that there is a 
universal repulsive core of the two-body force, but modified by the appearance of smeared DO-branes. And 
a short-distance attractive force would appear if () > 2, which consistently corresponds to the existence of 
unstable baryonic states in two-body system. The summary and conclusion are given in the final section. 


2 A brief review of Sakai-Sugimoto model in the D0-D4 background 

Here we begin to use D4’-brane to distinguish baryon vertex from those Nc D4-branes which are producing 
the background geometry. As we are going to derive the low energy effective theory for D4’-branes from 
Sakai-Sugimoto model in D0-D4 background, thus in this section, we first review the Sakai-Sugimoto model 
in D0-D4 background geometry briefly. Some of the results in this section are already presented in [251 [27]. 

2.1 D0-D4 background geometry 

By taking the near horizon limit, the solution of D4-branes with smeared DO charges in Type HA supergravity 
reads [25] 
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ds^ 


/ TJ \ 

( i? ) {U) dT^) 


(1) 


We have written this metric m in string frame and r is a periodic variable. With the near horizon limit we 
also have the formulas for dilaton, Ramond-Ramond (R-R) 4-form and the 2-form which is 

/ JJ \ A 1 

= ds dir; /2 = ^J^dU A dr; U = dC^ = Ber, (2) 

where we have 


A = 


{2Trlsy gsNp ^ 
W 4 V 4 ’ 


B = 


{2Trlsf Ncgs 

a;4 


Ho = l 


Ek. 

m' 


fX) = i-r 


(3) 


dVli, £4 and a ;4 = Stt^/S are the line element, the volume form and the volume of a unit 5"^. Ukk is the 
coordinate radius of the bottom of the bubble, and V 4 is the volume of D4-brane. Np and Nc are the numbers 
of DO- and D4-branes respectively. DO-branes are smeared in the x°, ...,x^ directions. The relations to the 
QCD variables are deformed as 


W = 




9s = 


Ukk = ImkkXIIHp {Ukk) ■ 


2Mkk ’ 

Here in order to keep the back reaction of DO-brane, we also required Np to be of order W as in 


(4) 


2.2 Embedded D8-branes and baryon vertex in D0-D4 background 

In the Sakai-Sugimoto model, flavors of the dual gauge field theory are introduced by Nf D 8 /D 8 -branes as 
probes embedded in the background geometry. Here we embed these Nf D 8 /D 8 -branes into the background 
described by as muss]. By taking the probe limit i.e. Nc,Np ^ Nf which makes the back reaction of 
the D 8 /D 8 -brane to the background negligible. The fermions created by open strings can be added to the 
fundamental representation of U{Nc) x 17 _r/l(IV/) which is treated as groups of chiral symmetry, while the 
gauge fields are in adjoint representation. 

We employ the viewpoint in Sakai-Sugimoto model, baryons have been provided as a D4’-brane wrapped on 
5^ which is called the baryon vertex imiiniini- Such D4’-branes have to attach the ends of W fundamental 
strings since the S'^ is supported by W units of a R-R flux in the supergravity solution. In this way, the 
baryon charge equals to N^ quark charge in the corresponding field theory. According to these arguments, 
baryons are also the D4’-branes wrapped on S'^ in our D0-D4 background geometry O- As we are going 
to study on the baryons, thus we will focus on the baryon vertex in the Sakai-Sugimoto model in D0-D4 
background in next sections. 


3 The matrix model from D0-D4/D8 system 

In this section, we derive the matrix quantum mechanics and we are going to use the action to study the 
baryon spectrum and the two-body interaction by an effective potential. As mentioned, the matrix model is 
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just the low energy effective action on the k D4’-branes embedded in ffavor Nf D8-branes in the geometry 
described by O- Here we first give our result i.e. the action of our matrix model in D0-D4/D8 system, and 
we leave the details for derivation in next two subsections. 


3.1 The action 

The action of our matrix model for baryons from D0-D4/D8 system is 


S = 


XNcMkk 

547r 


(1 + Tr I dt {DoX^f - ^ 


+ DoU!^DQUJia ~ 0 f ^ ~ 2 ^ ^ 


+ 


1 


4A2M4^ (1 + C)" 

+A^cTr / dtAo . 


i^y 


+ D ■ + D ■ 


(5) 


As we can see, ([S]) describes a quantum mechanical system with U (k) symmetry where k is the baryon 
number and the trace is taken over U (fc) adjoint representation. Here M = 1,2,3,4 and A = gy^Xc is 
the ’t Hooft coupling. The unique dimension-ful constant Mkk is defined as ®. D and Aq are auxiliary 
fields while the fields X^ and oj are the dynamical fields. C is defined as C = ^Qq/^kk- Note that all the 
fields are bosonic and this matrix model describes the A:-baryon system in D0-D4/D8 system according to 
the holographic principle. 

The matrix model ([S]) is a deformed matrix model of [53]. Without smeared DO-branes, i.e. setting 
C = 0, ® returns back to the matrix model in |24| . Therefore the matrix model ® is also a deformed 
ADHM matrix model as claimed in [24]. By integrating out of the auxiliary field D, it yields a potential of 
a commutator term as (Tr [X^X]) . So our matrix model looks also close to the BFSS Matrix theory [28] 
(which is understood as an effective description of M-theory) or the IKKT matrix model (30] if the last term 
and the mass term are absent. Note that the quadratic term of X'^ and w would be negative if (j > 2 which 
corresponds to a system with imaginary mass. It may be understood as the constraint for the stable states 
of baryons in D0-D4 system i.e. if C > 2, baryons may not be stable in this system. This viewpoint is exactly 
the same as m by using the instantoii views for baryons. Since we will use our matrix model to describe 
baryons in D0-D4/D8 system, thus only (j < 2 is considered here for stable baryons. 

As low energy effective theory, all the fields in our matrix model ® is similar to the matrix model of |24] . 
Therefore we employ the symmetry for our matrix model from |24| . And the representation of the fields is 
summarized in Table 1. The covariant derivative is defined as 


DqUJ 


DqX^ 


OqU! — iAouj 


doX^ -i[Ao,X^], 
Dquj = OqUJ iA^u). 


So the matrix model ® is a quantum mechanical system with the following symmetry 


U{k) X SU{Nf) X S'!? (3), 


( 6 ) 


5 






Fields 

index 

U{k) 

SU{Nf) 

SU{2) X SU{2) 

XM 

M = 1,2,3,4 

adj 

1 

(2,2) 

^ia 

a = 1,2\i = 1, 2...Nf 

adj 

fund 

(1,2) 

Aq 


adj 

1 

(1,1) 

Ds 

s = 1,2,3 

1 

1 

(1,3) 


Table 1: Fields in the matrix model 


Nf is the number of the flavors in QCD. The first symmetry U {k) is a local symmetry while the symmetry 
SU (Nf) X SO (3) is a global symmetry. SO (3) could be interpreted as the rotation symmetry in our space 
where the baryons live. If we embed the rotation symmetry as 

SO (3) C SO (4) ~ SU (2) X SU (2), 

then the symmetry of the action would be easier to understand. And the additional dimension corresponds 
to the holographic dimension. By the mass deformation, the SO (4) symmetry breaks down to the SO (3) 
symmetry. 


3.2 Derivation from Holography 

In D0-D4/D8 system, our concern is the D4’-branes wrapped on S'^ in the background (HJ, whose low energy 
effective theory is described by the matrix model ([S]). In AdS/CFT duality, such a D4’-brane is named as the 
“baryon vertex” |15) . which is responsible for creating or annihilating baryonic states in the dual held theory. 

Since the D4’-branes live inside the flavor branes, the action of the D4’-brane is related to the background 
geometry, the R-R flux and also affected by the presence of the flavor branes. So in the low energy effective 
theory of baryonic D4’-branes, the strings connecting the baryon vertices and the flavor branes provide the 
held Lo in the bi-fundamental representation. The location of the D4’-branes in the transverse directions 
is specified by the diagonal eigenvalues of the field in the adjoint representation. Thus these diagonal 
eigenvalues represent the locations of the baryons in our real 3-dimensional space. Due to the curved geometry 
and the flux, only the bosonic fields are kept here since the deformation between our matrix model and ADHM 
matrix model breaks the supersymmetry explicitly. So we will not care about the fermion part in our theory. 

As in [23], we assume that the S'^ dependence can be trivially reduced although the D4’-branes are inside 
the flavor branes and wrapped on S'^, that means a dimensional reduction with no dependence along S'^. 
Therefore the derived action is in the time dimension only i.e. depended on time only. So next, we will derive 
the low energy effective action of k brane vertices system by using standard technique in string theory. 

3.2.1 Derivation from DBI part 

Let us start from the action for a single D4’-brane which is 


Sdv 

= Sdbi + Scs , 


Sdbi 

= —Tda j (t'^e~'*’^J—det{GMN-\-2.'Ka'FMN), 


Scs 


(7) 
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Here we have used ^ to represent the coordinates on baryon vertex. The convenient coordinates are used 
here as in [111 [25] 


y = r cos 9 


Z = r sin0, 


9 = 


U 


1/2 

KK 




( 8 ) 


where the flavor D8-brane is located ai y = 0. We consider a stable D4’-brane situated at r = 0 wrapped on 
S"*, in this case the DBI action is 

Sdbi = f dtH^^^UR^ {R/Uf'^ V-Goo, (9) 

9s J 

where the induced metric is 


Goo — 


I) 


3/2 


H, 


1/2 




( 10 ) 


We have induced the metric on the worldvolume of D4’-branes to the low energy effective theory and the 
index of X® is / = 1, 2, 3, i.e. runs for 3-dimensioiial space, so we obtain 


Sdbi = — 


TDi'LOi 


5, 


— J dtHQ^^UR^^l - (doX^f - 


^Ukk (R 

9 U [u 


{doZf. 


( 11 ) 


Since we keep k baryons at short distance and only the low energy effective theory is considered here, so we 
just need to expand (HU for small Z and small X and define 


X4 

c 


2 

3 



Hk. 

UicK 


3/2 




yields a quadratic action 


( 12 ) 


Sdbi = 


XNcMkk 

27^? 




dt 


-l + -{doX^) --[^-2)Ml^{xy 


(13) 


The high orders of X, Z and their derivatives have been dropped off. 

The kinetic term and mass term for X in matrix action ([S]) is given by (fT51) . However this leaves the mass 
term for w. We have assumed that the mass term for uj is 


~ Y 2 ~ (14) 

This (ITTl) is a natural guess from a comparison with [21] • Since our matrix model is the deformation from the 
model in |24| , we simply set the mass term of w as 1 /4 times of the mass term of X'^ which is an assumption 
from |24| . 
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3.2.2 Commutator terms 


We will compute the commutator terms in the matrix action ([5]). First we expand the generic Dirac-Born- 
Infield (DBI) action of a Dp-brane Q to the quadratic order 

Sdbi ^ (2W)' ^Tdp j det GmnGmnGpqF^^^ . (15) 

The relevant relations from T-duality here is (27ra') Am = Gmm ■ We therefore have 


{2Tia'fGMNGpQF^PF^^ = 2G°°G,jDoX^DoX^ + 2G°°G,,DoZDoZ 

(zTra'j 

(zTra'j 

= -2(DoX^0'-^(l + C)"A2Mi^[X^,X^]\ (16) 

We have used the value of the metric at U = Ukk- The commutator term can be rewritten by an auxiliary 
field if we consider the following action 


S = c J dtTr 
Integrating out the field D, yields 

S = c J dtTi 

Comparing (fl^ with (fl^ . we have 


2 (27ra'f + D ■ 


\ f f 


S = 


XNcMkk ^^ 3/2 


547r 


(1 + C)^^^ dt 


S^TT^ 1 

4A2M4 ^ (1 + 0 ' 


(^oy + 3 ■ f^^X^^X^^ 


(17) 


(18) 


(19) 


3.2.3 Chern-Simons term 

Finally, let us consider the Chern-Simons part of the Dp-brane action. Since /4 = dG^ is given in ([2|), thus 
the Chern-Simons term in Q can be rewritten as 


Scs = ^ y ^3 A F 2 

— A _5_Tr / /ii/i2M3a/3(7 p 

~ 2Tr 2 ■ I 

= WTr f dtAo, (20) 


((21)1) is nothing but the last term in matrix action (fl^ . 
















4 Baryon spectrum 

In this section, we will use our matrix model (O to study baryon spectrum of D0-D4/D8 system for k = 1 
case, i.e the single baryon case. So first let us compute the Hamiltonian for a single baryon k = 1 with 
general Nf first. Since for fc = 1 the field X is a number, thus all commutators with X could be dropped. 
This leaves the terms of cj. 




XNcMkk 

547r 

XNcMkk 

547r 


(1 + Cf/^ 

(1 + Cf/^ 


r 

\ 36^2 

1 / 

- -A 

/ dt 

_4A2M^^ 

-f 

1 — 1 

Dj + D ■ 


/ 




A i j j* \ ( i 

4:LUiLU2 + [UJiLUi ) 


+ - 2a;ia;f 


( 21 ) 


The so-called ADHM potential is given by (I^Tl) . Basically we should solve it for the construction of instantons 
in the ADHM formalism, however it is equivalent to minimize the ADHM potential here. And since there is 
no dynamics for Aq in our matrix action ([5]), thus we can integrate it out and the terms including Aq become 




547r 


X dt 


I 


iAouj^douJai - idoUJiAoUJai + {A^f wfw, 


547r 


XMkk (1 -I- Q 


3/2 


At] 


( 22 ) 


So the equation of motion for Aq is 


By inserting (^51) to 


doCOai tOoUJ^ ^ai “t“ 2 AqCU^ ^ai T 


we obtain 


547r 


XMkk (1 -|- 


= 0 . 


(23) 




XNcMkk 

547r 


(i + cf'Y 


dt 


1 






547r 


AM, 


KK ( 1 +C) 


3/2 


(24) 


Then we use the definition of the momentum conjugate to the held u 


pa _ 


dS 


d (douji) 

XNcMkk 

547r 


(1 + C) 


3/2 


547r 


doUJi -- iujcdoujoc - idoUJcUioc + t —-^ ioj^ 

2q]uj^j\ XMkk + I 


.(25) 


Therefore we obtain the Hamiltonian 
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H = p^d^u^f-L 


XNcMkk 

547r 


(1 + C) 


3/2 


dowfdouf + ^ 


+ 


(1 + C)^ (4wiw^*w^a;f + (ujIujI*)^ + - 2ujIujT 


3®7r2 

1 f / 54:77 A ^ 


- Ym - Y - 7^3 + - doUjfuJaiY 

4uj°^u}at \\XMkkJ (1 + C) 


4.1 Single flavor 

For single flavor case, i.e. Nf = 1, we use the following ansatz for w 


Then the Hamiltonian (|26l) is 


— Pa^ ^d — Pd' 


(26) 


(27) 


H = 


XNcMkk 
5477 
1 / 


(1 + C) 


3/2 


277r y 

2p2 ^^XMkk) (1 + C) 


3 3 




■Op' 


(28) 


where we have used the definition 2p^ = pip* + P 2 P 2 - L®! us analyze the terms in (|28l) . In the soliton picture, 
the first term can be interpreted as a self-repulsion of the instanton which is induced by the Chern-Simons 
term with the path-integration of Aq. The second term comes from the mass term of the matrix model and 
the curved background geometry while the third term comes from the path-integration over auxiliary field D, 
which corresponds to the ADHM potential. Thus all the terms are physical and modified by the appearance 
of smeared DO-branes. 

It is easy to find that the Hamiltonian could be minimized by a nonzero p. In the large A limit, we can 
reduce the Hamiltonian (1^5)) to a linear formula by putting p = xX°‘M^Y. Here both x and a are constants. 
And similar as done in [53], every term in ((551) scales with large A limit. As a result, the second term in (1551) 
can be negligible, then the value for p to minimize this Hamiltonian is computed as 


and 


P = 2 “ 2 / 332 ^ 2 / 3 ;)^- 2 / 3 ^-^ ^^- 7/6 ^ ^- 29 ) 

Hrmn = 2-^ 77 -^/^X^/^N,MkK (1 + Cf'"" ■ (30) 


4.2 Two flavors 

In this subsection, we need to focus on a more realistic case which has two flavors. In order to eliminate 
the contribution from the ADHM potential, we first have to satisfy the ADHM constraints = 0, or 

equivalently, use the following ansatz 
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Nf Nf Nf Nj 

= XI \^&=l\^ = X hd=2r ■ (31) 

i —1 i —1 i —1 i —1 

The ADHM potential disappears if this condition is satisfied. And without loss of generality, this can be 
achieved by using the following choice 


wr' = 


■ , ,*=2 — 


(32) 


After including A'^-dependence, the Hamiltonian is 


H = 


XNcMkk 

547r 


+ l (l-^C)Mi^ M 


r 1 

( ) 

’ ' 1 ' 

_2p2 

(y XMkk / 

(1 + C)^ 3 


4x2 


This can be minimized at 


/ 1 \-l/4 

P = 2-1/437/4^1/2 a - 1/2 (1 - ic) (1 + 0“'^" . 


(33) 


(34) 


Then the minimum value of the Hamiltonian is 


^ _ NcMkk 

Ve 2^' 

In the soliton approach [IS1I2Z1, the value of p is the size of instantons. Obviously in our theory, we find that 
all the values are modihed by C which is related to the appearance of smeared DO-branes and (1351) would be 
totally complex if ^ > 2. 


4.3 Quantization 

Since we have obtained the results on the vacuum of the matrix model, we can then quantize the Hamiltonian 
for the A: = 1, Nf = 2 circumstance and the results should correspond to the baryon spectrum. We first 
rewrite the Hamiltonian used in [TB] since we are going to use the same tricks, so we have 


H = 


XNcMkk 

2 


547r 

5p2 

XMkk J 


(36) 


By comparing (1361) with (1331) . the energy spectrum in our system can be obtained easily by replacing the Q 
and u!p,ujz used in [TH] as 


5 / 1 \ 1/2 

Q -^Q (1 + C) j i^p 07 ^z) (1 + f 1 ~ i^-^) ■ 


(37) 


Then the mass formula for the baryon excitation is 
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M = Mo + (1 + C) 


(38) 


where Mq = ■ In order to fit the data from experimental values, we set ^ = 1.9303 with Nc = 3 for 

real QCD and set M^k = 945MeV to fit the mass of p meson. The baryon mass spectrum from our matrix 
model (1551) is listed in Table 2 which is more close to the experimental data . As a comparison, we also list 
experimental data in Table 3 and the superscripts ± represent the parity. The subscript * here is used to 
indicate that evidence for the existence of the baryonic states is poor. 

Such a baryon spectrum has already been obtained in |18j . However their original results are much larger 
than the experimental data if fitting the experiment data by setting Mkk = 945MeV (mass of p meson) 
with Nc = 3. And in |24| . the baryon spectrum would be still larger than the experimental data if fitting 
the experiment with the same value for Mkk and A).. Most likely, the reason is that Sakai-Sugimoto model 
describes the QCD with large limit by holography, but the real QCD is a theory with Nc = 3. Therefore, 
in our D0-D4/D8 system, we suggest to give an effective description for Nc = 3 QCD by adjusting the number 
density of DO-branes i.e. the parameter ( in our system. Note that our result (1551) does not make sense if 
^ > 2 since the mass spectrum would be totally imaginary. As mentioned, the stable baryonic state may not 
exist if C > 2 in D0-D4/D8 system, which would be quite different from the original Sakai-Sugimoto model. 


{np,n^) 

(0,0) 

(1,0) (0,1) 

(1,1) (2,0)/(0,2) 

(2,1)/(0,3) (1,2)/(3,0) 

TV (1 = 1) 

A (1 = 3) 

945+ 

1237+ 

1268+ 1268" 
1560+ 1560- 

1590- 1590+,1590+ 

1882- 1882+,1882+ 

1913-,1913- 1913+,1913+ 

2205-,2205- 2205+,2205+ 


Table 2: Baryon spectrum of mass from (1551) 


(»^p,^^^) 

(0,0) 

(1,0) (0,1) 

(1,1) (2,0)/(0,2) 

(2,1)/(0,3) (1,2)/(3,0) 

TV (1 = 1) 
A(1 = 3) 

940+ 

1232+ 

1440 1535- 

1600+ 1700- 

1655- 1710+,? 

1940- 1920+,? 

20907,? 2100+,? 

? ? ? ? 


Table 3: Experimental data of baryon mass for various baryonic states 


5 Two body baryon interaction 

In this section we will consider the interaction between two baryons for two-flavor case since it would be 
more realistic, thus TV/, k = 2. For the two-flavor case in the matrix model ([5]), integrating out of the U{k)- 
adjoint field Dab gives the vacuum configuration, which is just the ADHM constraints 

fA = 0 , (39) 

where the indices for baryon are A, B = l,2...fc. Since our matrix model is also a deformed ADHM matrix 
model, we will also use the ADHM data for our model as done in |24| 


Oil 


uJXi * — ^Pi 


2 

A=2 _ ttA=2 
ai ^ ai P2 


-X — T — + T Ym , 


(40) 
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where tm is the inter-baryon distance, is the SU (2) matrices which represents the moduli parameters of 
each baryon and 


Ym = 


-^i^Tr 

4(rL) 


CTMrNCTN 



(41) 


we have used ctm = 1) and ctm = {—iT, 1). Note that the ADHM data are just the solution for two Yang- 

Mills instantons which has been explicitly used in the approach of soliton [22]. So after the quantization, 
their degrees of freedom become the spin and the isospin of each baryon, which are nothing but the gauge 
rotations of the flavor gauge group. Here we use real unit vectors to write them as 


with {afY + {aff 
use 


= a f + iafY 


(42) 


1. Then we can obtain some useful expression listed as follows which we are going to 


tmYm = 0 , 

YmYm = -^^Tr 


8 {tm ) 


(c/2)i[7iy 


2 2 
P 1 P 2 

^{tmY 


1 ~ Wm^Ym) 


Tr 




= Tr 


=2a\,ai,. 


(43) 


Then, we shall compute the two-baryon interaction potential. First we write down the terms in the action 
(|5|) with ADHM data 


Tr (DoA“) = 2 + {Alf + 4Y^) + 4 (Ag) - ^AIaItmYm • (44) 

Note that the last term in vanishes once (H51) is added. Then, for the kinetic term of oj we have 


Tr (DowfDoWia) 


HpI + pD {Alf Y {Alf + {Alf + {Air 


+4pip2A°AjTr + 4 (p? _ ^2) ^0^, 


3 

0 • 


(45) 


We can minimize it with Aq = 0 since the component A§ appears only as (Aq)^. Thus the resulting baryon 
interaction potential H from the kinetic term plus Chern-Simons term can be evaluated by f dtV = —S'on-sheii, 
which is 


5 


on—shell 
kinetic 


XNcMkk 

547r 

XNcMkk 

547r 


(1 + C)^/" Tr J dt [{DoX^) + DoojfDoUj^c 

(1 + I dt 2 {Al)\l, + 8 (A^) 


+ NcTr / dtAo 


+ 2 (pj -I- P 2 ) (^(Ag) + (Aj) -I- (Ag) ^ 

(([/■)' U^) + 4 W - pS} AlAl + (1 + 0-“''“ Ap 


(46) 
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Since the action (1461) does not depend on the derivative of Aq, we can integrate them out straightforwardly, 
therefore the resulting baryon potential is 


V = 


277rNc 


(pf + pI + rli) {upIpI -A{pI+ pI) r 


XMkk (1 + () 3/2 pIpI 


u (pf - (2 + u) pIpI + p\) + 5u {pI + pI) rh - 16 {rhY 


(47) 


where we have defined u = 


(Tr (C/i)^ 




— 4. And there is another part of baryon interaction potential, 


which comes from the mass term of in action ([51) in addition to this potential (|47p . It can be evaluated as 


(1 + C)"/" I (1 - K) ^KKTr (x^f 

= + • (48) 

This term (H51) does not contribute to the two-baryon interaction since the term of is the mass term in the 
single-baryon Hamiltonian. In fact, the off-diagonal elements of Y 4 contribute to the interaction between the 
baryons. Then using (HIT) we have 


T4 = - 


P 1 P 2 

M 




where i = 1,2,3, so we can rewrite the potential energy 


as 


(49) 


XNcMkk 

1627r 


(1 + C) 


3/2 




2 

KK 


pIpI 

{rltf 


r,:Tr 




Therefore, putting all together, we obtain the two-body interaction Hamiltonian of baryons which is 


V = 


277rNr 


(pf +pI+ rh) (upfpl -4 (pI + pI) r|f) 


\Mkk (i + ^)3/2 p\p\ 


I ^^cMkK ^ ^^3/2 


1627r 

27TTNr 


U (pi - (2 + u) pIpI + pY) + 5m (pI + pI) rlj - 16 (rljY 

1 \ 2 


1 - oC ) 


PipI 

(rl,f 


r,;Tr 


■.^i^(l)t[/(2) 


4AMkk (1 + ^)3/2 \p{ p2 


1 1 
3/2 1 72 + 72 


(50) 


(51) 


If we assume that the size p of the baryons or nucleons is much small, then we can expand (1511) for tm p 
and obtain a leading term which is 


V = 


277rNr 


1 


1 


64XMkk (1 -h C)^^^ I 


32 -|- 6m -I- (5u -I-16) 1 

\Pi 


2 “T 9 


pfyj 






i+uAn[/(2) 


1 \ 2 


(52) 


again by choosing pi = P 2 = P, we have a formula from (1521) 


14 



































V = 


277rNr 


AXMkk (1 + rlj 


^Tr 




ir%) 


(r.Ti 


(53) 


Then, let us compute the vacuum expectation value by using the potential (1531) . We use ^/aj Xii ''^pi 
to label the states of the two baryons, with A = 1,2 representing the two baryons. The excited baryon states 
are labeled by Up and nz which are the quantum numbers. And /, J is the isospin (spin) of the baryon state 

respectively. For nucleons we have ^ / = J = 5 ^ the explicit spin/isospin wave functions reads 


\pt> \pi> 

|nt> \ni> 


(54) 


The baryon potential is evaluated as 


16 


h,Jl,l2,J2 


= 6 ^^ + y Jl • I2 + Ji J\ - ■ J2) 


Tr Tr 

By the standard definition of S '12 = 12 • rj — dj) • J 2 with r = r/ \r\, we have 


(55) 


{r) = TT 


I+8 (/../,) (/../,) 

q“'(rl = 2»(/, ./j) 


N, 


1 1 


do) 


XMkk (1 + (C)^/^ r"^ ’ 

1 


XMkk (1 + ' 

For the leading order which comes from the second term in (l53l) . we hav^l 


(56) 


yW(r) = 


yW(r-0 = 


1 16 


81 2187 

8 


(/i • h) (Ji • J 2 ) 


3 / 2 /, i;\\N,Mij,p^ 




2187 


(/../.) ( 1 +<r(i-if) 


Thus our total potential up to leading order is 


9 ’ 


(57) 


Vc{f) = (r) + 1/W (rl, 

Vrif) = (r) + 1/W ^ 

Vnuclear = Vq {r) + S12VT (F) . 


^Here the exact result is 




■J , we simply write it as 


(58) 
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This is the short-distance two-body force for baryons in D0-D4/D8 system obtained from our matrix 
model. With ^ < 2 we also find there is a repulsive core of baryons or nucleons in this system. As in [2^ 124) . 
the repulsive potential scales as r~^ which was treated as a property peculiar to the holographic model . If 
comparing ((53l) ((56l) ((57l) with the results in [24| , we can see that the nuclear force in |24] is just our results 
with C = 0 i.e. no smeared DO-branes. With the appearance of smeared DO-branes, the zeroth-order result is 
depressed by the number of DO while the leading order result increases. As mentioned in previous sections, 
our model for baryons is consistent with C < 2 and this is the condition for the existence of stable baryons 
in our theory. However an attractive force appears from (1571) if C > 2 which may give an unstable state for 
baryons in two-body system by this model. 


6 Summary and discussion 

We have studied the holographic baryons by the matrix model in the D0-D4/D8 system. We start with 
considering the baryon vertex inside the flavor branes i.e. D4’-branes wrapped on S'^ which is embedded in 
the D0-D4 background on large Nc limit. And we use the standard technique in string theory to derive our 
matrix model. By using T-duality and dimensional reduction, we obtain our matrix model ([5]) with U (fc) 
symmetry, which could also be able to describe multi-baryons. In order to describe stable baryonic states, 
we find the value of C is restricted to ^ < 2. This is exactly the same as the results in m- However in m, 
baryons are described by using BPST configurations which is similar as [TB]. In our paper, as a difference, 
we start from the baryon vertex, but we come to the same conclusion i.e. stable baryons exist in D0-D4/D8 
system only if C < 2. So it is worthy to believe such a result is unique for this system. 

With this matrix model, we also determine the holographic size and the spectrum of baryon for the case 
with k = 1 and Nf = 1,2. We find the spectrum obtained from our model could be more close to the 
experimental data by choosing suitable value of Thus we suggest that it is an effective description of QCD 
for low energy baryonic states. Again, we have seen that in the two flavor case the spectrum of baryon’s 
mass would be totally imaginary if C > 2. It turns out that baryons in this system are not stable if ^ > 2 as 
mentioned. In two flavor case, our spectrum of baryon’s mass (1551) could recover the results in [53] if setting 
C = 0 i.e. no smeared DO-branes. However, we can not recover the results (l36l) from [T8|, since the first term 
in ((33l) is deformed by a factor 4/5, even if setting C = 0- This puzzle was found in [24] first which makes our 
spectrum (I38L from the matrix model of D0-D4/D8 system, a bit different from the results by the instanton 
viewpoint in m- We are less clear about this and would take further study about this in the future work. 

And we also compute the two-body short-distance effective potential for baryons in D0-D4/D8 system 
i.e. k = 2 case. It also exhibits a repulsive core and a tensor force as [221124| . which has been well-known in 
nuclear physics, but as a difference, it is modified by the appearance of smeared DO-branes. As it is known 
the nucleons interact each others by interchanging mesons and the effective potential can be derived from 
the Yukawa coupling. And on the other hand, in [251 ED, the spectrum of mesons, such as p mesons, of 
Sakai-Sugimoto model in D0-D4 background are modified by the appearance of smeared DO-branes. So in 
our paper, (1551) (1571) can be interpreted as follows: since the mass of mesons are modified, thus the effective 
potential from the interaction are also modified by the appearance of smeared DO branes. This two-body 
interaction has not been calculated by using the instanton viewpoints in D0-D4/D8 system in m , so our 
work would be advancing on this front. Our effective two-body potential (1551) (157)) is obtained by expanding 
((5T1) respected to p/r, however it shows that if C>2, the leading order potential would be negative i.e. it is an 
attractive force. Furthermore, if C and A are large enough, the zero order potential would be depressed very 
much while the leading order potential becomes a rapid negative increasing, which makes the total potential 
negative. A negative two-body short-distance effective potential implies there is an attractive force between 
baryons or nucleons at short distance, so the system would not be stable. From the effective two-body force, it 
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also and consistently turns out that is restricted even in two-body system if there would be stable baryonic 
states. 

Finally, we would like to give some more comments about this work. Basically, our analysis is done in large 
Nc limit as in most analysis in gauge-gravity duality, though we have tried to give an effective description 
of finite Nc theory in this paper. The interaction between baryons would be important for studying on the 
phase structure of strong-coupling QCD by holography, such as imiini. So a holographic model corresponded 
to real QCD would be significant, however unfortunately, currently out of reach. Besides, we can see that 
our two-body effective potential for baryons makes sense at short distance only, and furthermore we still 
do not know how to introduce a realistic attractive force consistently or describe the long-distance force by 
considering the baryon vertex directly. These problems may be solved by calculating the interaction between 
the D4’-branes in curved spacetime, however it would be a huge challenge in string theory. So it should be 
understood as that there is still a long way to the realistic and strong-coupling QCD from holography. 
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